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Abstract 

We define a higher spin ahernating sign matrix to be an integer-entry square 
matrix in which, for a nonnegative integer r, all complete row and column 
sums are r, and all partial row and column sums extending from each end 
of the row or column are nonnegative. Such matrices correspond to config- 
urations of spin r/2 statistical mechanical vertex models with domain- wall 
boundary conditions. The case r = 1 gives standard alternating sign matri- 
ces, while the case in which all matrix entries are nonnegative gives semimagic 
squares. We show that the higher spin alternating sign matrices of size n are 
the integer points of the r-th dilate of an integral convex polytope of dimen- 
sion (n— 1)^ whose vertices are the standard alternating sign matrices of size n. 
It then follows that, for fixed n, these matrices are enumerated by an Ehrhart 
polynomial in r. 
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1. Introduction 



Alternating sign matrices are mathematical objects with intriguing combinatorial 
properties and important connections to mathematical physics, and the primary 
aim of this paper is to introduce natural generalizations of these matrices which also 
seem to display interesting such properties and connections. 

Alternating sign matrices were first defined in [IH] , and the significance of their con- 
nection with mathematical physics first became apparent in [16] , in which a determi- 
nant formula for the partition function of an integrable statistical mechanical model, 
and a simple correspondence between configurations of that model and alternating 
sign matrices, were used to prove the validity of a previously-conjectured enumera- 
tion formula. For reviews of this and related areas, see for example [Ul [ISl EH [7T] . 
Such connections with integrable statistical mechanical models have since been used 
extensively to derive formulae for further cases of refined, weighted or symmetry- 
class enumeration of alternating sign matrices, as done for example in [211 HTf [55 | [70] . 

The statistical mechanical model used in all of these cases is the six-vertex model 
(with certain boundary conditions), which is intrinsically related to the spin 1/2, 
or two dimensional, irreducible representation of the Lie algebra s/(2,C). For a 
review of this area, see for example [38]. In this paper, we consider configurations 
of statistical mechanical vertex models (again with certain boundary conditions) 
related to the spin r/2 representation of sl{2, C), for all nonnegative integers r, these 
being in simple correspondence with matrices which we term higher spin alternating 
sign matrices. Determinant formulae for the partition functions of these models 
have already been obtained in [T7j, thus for example answering Question 22 of [17] 
on whether such formulae exist. 

Although we were originally motivated to consider higher spin alternating sign ma- 
trices through this connection with statistical mechanical lattice models, these ma- 
trices are natural generalizations of standard alternating sign matrices in their own 
right, and appear to have important combinatorial properties. Furthermore, they 
generalize not only standard alternating sign matrices, but also other much-studied 
combinatorial objects, namely semimagic squares. 

Semimagic squares are simply nonnegative integer-entry square matrices in which 
all complete row and column sums are equal. They are thus the integer points 
of the integer dilates of the convex polytope of nonnegative real-entry, fixed-size 
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square matrices in which all complete row and column sums are 1, a fact which 
leads to enumerative results for the case of fixed size. For reviews of this area, see 
for example [7, Ch. 6] or [62l Sec. 4.6]. In this paper, we introduce an analogous 
convex polytope, which was independently defined and studied in [M], and for which 
the integer points of the integer dilates are the higher spin alternating sign matrices 
of fixed size. 

We define higher spin alternating sign matrices in Section 2, after which this pa- 
per then divides into two essentially independent parts: Sections 3, 4 and 5, and 
Sections 6, 7 and 8. In Sections 3, 4 and 5, we define and discuss various combi- 
natorial objects which are in bijection with higher spin alternating sign matrices, 
and which generalize previously-studied objects which are in bijection with standard 
alternating sign matrices. In Sections 6, 7 and 8, we define and study the convex 
polytope which is related to higher spin alternating sign matrices, and we obtain 
certain enumerative formulae for the case of fixed size. We then end the paper in 
Section 9 with a discussion of possible further research. 

Finally in this introduction, we note that standard alternating sign matrices are 
related to many further fascinating results and conjectures in combinatorics and 
mathematical physics beyond those already mentioned or directly relevant to this 
paper. For example, in combinatorics it is known that the numbers of standard 
alternating sign matrices, descending plane partitions, and totally symmetric self- 
complementary plane partitions of certain sizes are all equal, but no bijective proofs 
of these equalities have yet been found. Moreover, further equalities between the 
cardinalities of certain subsets of these three objects have been conjectured, some 
over two decades ago, and many of these remain unproved. See for example [3l HI 
ESI ESI Hni SD Eni EI] • Meanwhile, in mathematical physics, extensive work has been 
done recently on so-called Razumov-Stroganov-type results and conjectures. These 
give surprising equalities between numbers of certain alternating sign matrices or 
plane partitions, and entries of eigenvectors related to certain statistical mechanical 
models. See for example [23l [23] and references therein. 

Notation. Throughout this paper, P denotes the set of positive integers, N denotes 
the set of nonnegative integers, [m, n] denotes the set {m, m -|- 1, . . . , ra} for any 
m, n e Z, with [m, = for n < m, and [n] denotes the set [l,n] for any n G Z. 
The notation M(o,i) and M[o,i] will be used for the open and closed intervals of real 
numbers between and 1. For a finite set T, |T| denotes the cardinality of T. 
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2. Higher Spin Alternating Sign Matrices 



In this section, we define higher spin alternating sign matrices, describe some of 
their basic properties, give an enumeration table, and introduce an example. 

For n G P and r G N, let the set of higher spin alternating sign matrices of size n 
with line sum r be 
ASM(n,r) : = 

/An . . 



A-- 



Aln\ 



G 



\Ani ■ ■ ■ Ann) 



• E;LiA,,' = T.7'=iAvj = r for all 1,3 G [n]\ 

• Ej'=i Aj' > for all i G [n], j G [n-1] 

• T.}=j Aij' > for all i e [n], j e [2, n] 

• Td'=iAi'j > for all i e [n-1], j e [n] 

• Y:7'=iAi'j > for all i G [2,n], j G [n] 



>. 



In other words, ASM(n, r) is the set of nxn integer-entry matrices for which all 
complete row and column sums are r, and all partial row and column sums extend- 
ing from each end of the row or column are nonnegative. As will be explained in 
Section 3, a line sum of r corresponds to a spin of r/2. The set ASM(n, r) can also 
be written as 
ASM(n,r) = 

E"'=i Aj, = E?=i A,j = r for all t,j G [n] ' 
[A= ■■ ■■ G Z"^" • < Ej'=i Aij, < r for all i G [n], j G [n-1] >. 

• < Ei'=i Ai'j < r for all i G [n-1], j G [n] ^ 





■ Ain\ 








G Z"^" 


\Anl ■ 


■ Ann) 





(2) 



It follows that each entry of any matrix of ASM(n, r) is between — r and r, and that 
if the entry is in the first or last row or column, then it is between and r. 



A running example will be the matrix 

/O 1 1 0\ 

1-10 2 
A = 1 1-2 2 

10 10 

\0 1 1 0/ 



G ASM(5,2). 



(3) 



Defining 



SMS(n,r) := {A G ASM(n,r) | A^ > for each i,j G [n]}. 



(4) 
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it can be seen that this is the set of semimagic squares of size n with hne sum r, 
i.e., nonnegative integer-entry nxn matrices in which all complete row and column 
sums are r. For example, SMS(n, 1) is the set of nxn permutation matrices, so that 

|SMS(n,l)| = n\ (5) 

Early studies of semimagic squares appear in [21 HH]. For further information and 
references, see for example [TJ Ch. 6], [31], [60], [61], [621 Sec. 4.6] and [63l Sec. 5.5]. 

It can also be seen that ASM(n, 1) is the set of standard alternating sign matrices of 
size n, i.e., nxn matrices in which each entry is 0, 1 or —1, each row and column con- 
tains at least one nonzero entry, and along each row and column the nonzero entries 
alternate in sign, starting and finishing with a 1. Standard alternating sign matrices 
were first defined and studied in [IHl [50]. For further information, connections to 
related subjects, and references see for example [151 HSl [211 ISll 1561 [71]. 

We refer to ASM(n, r) as a set of 'higher spin alternating sign matrices' for any 
n G P and r G N, although we realize that this could be slightly misleading since the 
'alternating sign' property applies only to the standard case r = 1, and the spin r/2 
is only 'higher' for cases with r > 2. Nevertheless, we still feel that this is the most 
natural choice of terminology. 

Some cardinalities of ASM(n, r), many of them computer-generated, are shown in 
Table [H 





r = 


1 


2 


3 


4 


n=l 


1 


1 


1 


1 


1 


2 


1 


2 


3 


4 


5 


3 


1 


7 


26 


70 


155 


4 


1 


42 


628 


5102 


28005 


5 


1 


429 


41784 


1507128 


28226084 


6 


1 


7436 


7517457 


1749710096 


152363972022 



Table 1: |ASM(n,r)| for n e [6], r e [0,4]. 



Apart from the trivial formulae |ASM(0,n)| = 1 (since ASM(0,n) contains only the 
nxn zero matrix), |ASM(l,r)| = 1 (since ASM(l,r) = {(r)}), and |ASM(2,r)| = 

i G [0, r] [ = SMS (2, r)), the only previously- 



r+1 (since ASM(2, r) = {(^1, 



r — t 

i 



known formula for a special case of |ASM(n,r)| is 



|ASM(,U)I = nL_i:. ,6) 

for standard alternating sign matrices with any n G P. This formula was conjectured 
in |l9l[50], and eventually proved, using different methods, in [69] and [16]. It has 
also been proved using a further method in [34j , and, using a method related to that 
of [16], in [21]. 



3. Edge Matrix Pairs and Higher Spin Vertex Model Configurations 

In this section, we show that there is a simple bijection between higher spin alter- 
nating sign matrices and configurations of higher spin statistical mechanical vertex 
models with domain-wall boundary conditions, and we discuss some properties of 
these vertex models. 

For n G P and r G N, define the set of edge matrix pairs as 
EM(ra,r) := 



/ / HiQ . . . Hin\ 



(H, V) : 



fVoi . . . VonW 



^ l^g ^jnx(nrf-l) ^ ^j(rH-l)xn 



\\-f^nO • • • HnnJ • • • ynnj / 

Hio = Voj = 0, = Vnj = r, Hij_i + Vij = Vi^ij + Hij, for all i,j G [n] 



(7) 



We shall refer to if as a horizontal edge matrix and V as a vertical edge matrix. It 
can be checked that there is a bijection between ASM(n, r) and EM(n, r) in which 
the edge matrix pair {H, V) which corresponds to the higher spin alternating sign 
matrix A is given by 

i 

Hij = ^ Aijf , for each i E [n], j E [0, n] 



Vij = ^ Ai'j , for each i G [0, n], j G [ 



n\ 



i'=i 



and inversely. 



Aij = Hij - Hij_i = Vij - Vi^ij , for each i,j e [n] 



(9) 
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Thus, H is the column sum matrix and V is the row sum matrix of A. The cor- 
respondence between standard alternating sign matrices and edge matrix pairs was 
first identified in [58] . 

It can be seen that for each {H, V) G EM(n, r) and i,j G [0, n], I]"=i Hi'j = jr and 
J2]'=i Vij' = ir, so that 

n n 

J2 H,, = J2 = n{n + l)r/2. (10) 

i,j=l i,j=l 



The edge matrix pair which corresponds to the running example ([3]) is 



iH,V) 



/O 1 2 

10 

12 

111 

yo 1 1 



2\ 

2 
2 
2 

2/ 



/O o\\ 
110 
10 12 
112 2 
2 12 12 

V2 2 2 2 2/ / 



A configuration of a spin r/2 statistical mechanical vertex model on an nxn square 
with domain-wall boundary conditions is the assignment, for any {H, V) G EM(n, r), 
of the horizontal edge matrix entry Hij to the horizontal edge between lattice points 
{i,j) and (i, j+1), for each i G [ri], j G [0, n], and the vertical edge matrix entry Vij to 
the vertical edge between lattice points (z, j) and j), for each i G [0, n], j G [ra]. 
Throughout this paper, we use the conventions that the rows and columns of the 
lattice are numbered in increasing order from top to bottom, and from left to right, 
and that denotes the point in row i and column j, i.e., we use matrix-type 

labeling of lattice points. The assignment of edge matrix entries to lattice edges 
is shown diagrammatically in Figure [H and the vertex model configuration for the 
example of (fTTj) is shown in Figure [2l The term domain-wall boundary conditions 
refers to the assignment of to each edge on the left and upper boundaries of the 
square, and of r to each edge on the lower and right boundaries of the square, i.e., 
to the conditions Hio = Voj = and ifj„ = Vnj = r of ([7]). The correspondence 
between standard alternating sign matrices and configurations of a vertex model 
with domain- wall boundary conditions was first identified in |32j . 

We note that in depicting vertex model configurations, it is often standard for certain 
numbers of directed arrows, rather than integers in [0,r], to be assigned to lattice 
edges. For example, for the case r = 1, a configuration could be depicted by assigning 
a leftward or rightward arrow to the horizontal edge from {i,j) to (i, for Hij = 
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Figure 1: Assignment of edge matrix entries to lattice edges. 
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Figure 2: Vertex model configuration for the running example. 



or Hij = 1 respectively, and assigning a downward or upward arrow to the vertical 
edge between and + for Vij = or = 1 respectively. The condition 
ifi j_i + Vij = Vi_i J- + Hij of ([7]) then corresponds to arrow conservation at each 
lattice point (i.e., that the numbers of arrows into and out of each point are equal), 
while the domain-wall boundary conditions correspond to the fact that all arrows 
on the horizontal or vertical boundaries of the square point inwards or outwards 
respectively. 

It is also convenient to define the set of vertex types, for a spin r/2 statistical 
mechanical model, as 

V(r) := {(/i, V, h', v') E [0, r]^ | h + v = h' + v'}. (12) 



>h', and it can be seen that for the 



A vertex type {h, v, h', v') is depicted as h*- 

vertex model configuration associated with (if, V) G EM(n, r), the lattice point («, j) 
is associated with the vertex type {Hi j_i, Vij, H^j, Vi_ij) G V(r), for each i,i G [n]. 
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(16) 
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(17) 



1 

(18) 
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(19) 



Figure 3: The 19 vertex types of V(2). 



The vertex types of V(2) are shown in Figure [3], where (1)-(19) will be used as 
labels. The vertex types of V(l) are (l)-(5) and (10) of Figure [3l 

For any r G N, V(r) can be expressed as the disjoint unions 

2r 

y{r) = [J \{h, s — h,h' , s — h') /;,,/;,' G [max(0, s—r), min(r, 



s=0 



{{h, V, h', h+v~h') I h, v, h' G [0, r], h < h' < v} U 
{{h, f , h+v-v', v') I h, V, v' G [0, r], v < v' < h} U 
{{h, h'+v'-h, h', v') I h, h', v' G [0, r], h' < h < v'} U 
{{h'+v'-v,v,h',v') I v,h',v' G [0,r], v' < v < h'} , 



so that 



(13) 



(r+l)(2r2+4r+3)/3. (14) 



It can be seen, using (jlj) and (Q, that a spin r/2 vertex model configuration cor- 
responds to a semimagic square with line sum r if and only if each of its vertex 
types is in Vs(r) := {{h,v, h' ,v') G V(r) \ h < h' (and v' < v)}. For example, Vs(l) 
consists of (l)-(3), (5) and (10) of Figure El and Vs(2) consists of (l)-(3), (5)-(8), 
(10), (11), (14)-(16), (18) and (19) of Figure! 
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By imposing the condition h < h' on the two disjoint unions of f|T3l) . which in the 
second case leaves just the first and fourth sets, it follows that |Vs(r)| = J^g^i^"^ = 
r + 2\ /r + 3^ 



(r + l)(r + 2)(2r + 3)/6. 

For a spin r/2 statistical mechanical vertex model, a Boltzmann weight 

W{r,x,h,v,h\v') eC (15) 

is defined for each {h,v, h' ,v') G V(r). Here, x is a complex variable, often called 
the spectral parameter. 

For such a model on an n by n square with domain-wall boundary conditions, and 
an nxn matrix z with entries Zij G C for i,j G [n], the partition function is 

n 

Z{n,r,z) := E n W{r, z,„ K,,.,,V,„ H,„V^,j) . (16) 

{H,V)eEM{n,r) i,j=l 

Values of Z{n,r, z) therefore give certain weighted enumerations of the higher spin 
alternating sign matrices of ASM(n, r). It follows that if there exists G C such 
that 

W{r,u^,h,v,h',v') = 1 for each (/i, t>, /i', t>') G V(r), (17) 

then 

Z{n,r,z)\^^^i,,^^=^A = |ASM(ra,r)|, (18) 
and that if there exists u^. G C such that 

W{r, <, h, V, h', v') = |J' ^ ^, for each {h, v, h', v') G V(r), (19) 

then 

^(n,r,^)Uach.„=«s = |SMS(n,r)|. (20) 



The Boltzmann weights f|T5|) are usually assumed to satisfy the Yang-Baxter equation 
and certain other properties. See for example pi Ch. 8 & 9] and [381 Ch. 1 & 2]. Such 
a model is then known as integrable, and is related to the spin r/2 representation, 
i.e., the irreducible representation with highest weight r and dimension r+l, of the 
simple Lie algebra s/(2,C), or its affine counterpart. See for example [351 EHl 138] . 
Each value i G [0,r], as taken by the edge matrix entries, can thus be associated 
with an sl{2, C) weight 2i—r. In physics contexts, it is also natural to associate each 
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i G [0, r] with a spin value i—r/2. The model with r = 1 is known as the six-vertex 
or square ice model, and is related to the XXZ spin chain and the defining spin 1/2 
representation of sZ(2,C). Furthermore, the Boltzmann weights for each case with 
r > 1 can be obtained from the Boltzmann weights for the r = 1 case using a 
procedure known as fusion. See for example [15] . Boltzmann weights for the r = 2 
case are also obtained more directly in [391 EH EH] ■ 

For the Boltzmann weights of these models, and for any x = {xi, . . . ,Xn), y = 
(j/i, . . . , Un) G C" with each having distinct entries, it can be shown that 

Z(ra,r,2;)|each2i,=xi^, = F{n,r,x,y) det M(n, r, x, , (21) 

where M{n,r,x,y) is an nr x nr matrix with entries M{n,r,x,y)(^i^k),{j,i) = (f){k — l, 
Xi—yj) for each {i,k), (j, /) G [n] x [r], and F and are relatively simple, explicitly- 
known functions. This determinant formula for the partition function with Zij = 
Xi—yj is proved for the case r = 1 in [ISl |33], using results of [S], and for each case 
with r > 1 in ^7\, using the r = 1 result and the fusion procedure. The formula for 
the r = 1 case is also proved in [TT], using a method different from that of [121 
while that for each case of r > 1 was obtained independently of [17], but using a 
similar fusion method, in [8]. 

If any entries of x, or any entries of y, are equal, then F{n, r, x, y) has a singularity, 
and det M(n, r, X, ?/) = 0. However, by taking an appropriate limit as the entries 
become equal, as done in [43] for r = 1 and [I^ for r > 1, a valid alternative formula 
involving the determinant of an nr x nr matrix whose entries are derivatives of the 
function can be obtained. For the completely homogeneous case in which all entries 
of X are equal, and all entries of y are equal, with a difference u between the entries 
of X and this matrix has entries (pik — l, u) for each (z. A;), (j, /) G [n] x [r]. 

For the case r = 1, there exists such that the Boltzmann weights which lead to 
a determinant formula fl^ satisfy ffTTI) . so that flTSl) can be applied. This is done 
in [16] and [21] in order to prove ([6]). In ^46j, a choice of x and y which depend 
on a parameter e is used, in which x and y each have distinct entries for e 7^ 0, 
and Xi — yj = for e = and each i,j G [n]. The formula (12T|) is then applied 
with e 7^ 0, the resulting determinant is evaluated as a product form, and finally 
the limit e — > is taken, giving the RHS of ([6]). In [21], a determinant formula for 
the completely homogeneous case is applied at the outset, and the relation between 
Hankel determinants and orthogonal polynomials, together with known properties of 
the Continuous Hahn orthogonal polynomials, are then used to evaluate the resulting 
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determinant, giving the RHS of IQ. 

For the cases with r > 1, if there exist values and sucli that Boltzmann weights 
which lead to a determinant formula (12T1) satisfy (|T71) and (fT9l) . then methods similar 
to those used in the r = 1 case could be applied in an attempt to obtain formulae 
for |ASM(n, r)| and |SMS(n,r)| for fixed r and variable n. However, unfortunately 
it seems that such and might not exist. 

4. Lattice Paths 

In this section, we show that there is also a bijection between higher spin alternating 
sign matrices and certain sets of lattice paths. 

For G P and r G N, let LP(?7,, r) be the set of all sets P of nr directed lattice paths 
such that 

• For each i & [n], P contains r paths which begin by passing from 
(n+l,i) to {n,i) and end by passing from {i,n) to {i,n + l). 

• Each step of each path of P is either ( — 1, 0) or (0, 1). 

• Different paths of P do not cross. 

• No more than r paths of P pass along any edge of the lattice. 

It can be checked that there is a bijection between EM(n, r) (and hence ASM(?t,, r)) 
and LP(n, r) in which the edge matrix pair {H,V) which corresponds to the path 
set P is given simply by 

Hij = number of paths of P which pass from 
to + for each i G [n], j G [0,n] 

(22) 

Vij = number of paths of P which pass from (i + 1, j) 
to (i, j), for each i G [0,n], j G [n]. 

For the inverse mapping from [H, V) to P, fl22l) is used to assign appropriate numbers 
of path segments to the horizontal and vertical edges of the lattice, and at each 
(z, j) G [n] X [n], the Hij_i+Vij = Vi^ij+Hij segments on the four neighboring edges 
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are linked without crossing through according to the rules that 

• If Hij = Vij (and -ffij-i = ^-ij), then i/ij-i paths pass from — 
to (i — 1, j), and Hij paths pass from + j) to (i, j + 1). 

• If Hij > Vij (and Hi j_i > Vj^ij), then Vj^ij paths pass from (i,j — 1) 

to (i - 1 , j) , - l/jj = Hij_i - Vi^ij paths pass from ^^3^ 
to (i, j + 1), and V^^ paths pass from (i + 1, j) to (i, j + 1). 

• If Vij > Hij (and l^-ij > Hij^i), then paths pass from (i, j — 1) 
to (i — l,j), Vij — Hij = Vi^ij — Hi j^i paths pass from (i + l,j) 

to (i — 1, j), and Hij paths pass from (i + 1, j) to (i, j + 1). 

The three cases of (1231) are shown diagrammatically in Figure HJ the path config- 
urations which correspond to the vertex types of V(2) from Figure [3] are shown in 
Figure [5l and the path set of LP(5,2) which corresponds to the running example 
of ([3]), (ITT!) and Figure [2] is shown in Figure [61 In order to assist in their visualiza- 
tion, some of the path segments in these diagrams have been shifted slightly away 
from the lattice edges on which they actually lie. Also, as indicated in the previous 
section, we are using matrix-type labeling of lattice points. 



(i,j-l)G- 



[i- 

C 


,j) 

p 


J. 






'( 




■ "ij 




Hij -- 





3(i,ifl) {i,J-l)o- 



(i.i) 



Hij > Vij 



(i,j-l)G- 



Vij > Hij 



Figure 4: Path configurations through vertex {i,j) for the cases of (l23l) . 



The case LP(n, 1) of path sets for standard alternating sign matrices is studied in 
detail in [9] as a particular case of osculating paths which start and end at fixed 
points on the lower and right boundaries of a rectangle. The correspondence between 
standard alternating sign matrices and such osculating paths is also considered in 
Sec. 5], [H Sec. 2], [30l Sec. 9] and [65, Sec. IV]. 
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Figure 5: Path configurations for tlie 19 vertex types of V(2). 
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Figure 6: Set of lattice paths for the running example. 

5. Further Representations of Higher Spin Alternating Sign Matrices 

In this section, we describe three further combinatorial objects which are in bijection 
with higher spin alternating sign matrices: corner sum matrices, monotone triangles 
and complementary edge matrix pairs. These provide generalizations of previously- 
studied combinatorial objects which are in bijection with standard alternating sign 
matrices. 
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For n G P and r G N, let the set of corner sum matrices be 



g |^(n+l)x(nfl) 



CSM(n,r) := (C- 

















\CnO ■ 


■ Cnn J 



(24) 



Cok = Cko = 0, Ckn = Cnk = kr, for all k G [n] 

< < r, <Cij-Q^ij <r, for alH, j G [n] 



It can be checked that there is a bijection between ASM(n,r) and CSM(n, r) in 
which the corner sum matrix C which corresponds to the higher spin alternating 
sign matrix A is given by 



t 3 



and inversely, 



A, 



i'=lj'=l 



Cij - Cij^i - Ci^ij + C^ij-i , for each i,j G [n]. 



(25) 



(26) 



Combining the bijections (151151) between EM(n, r) and ASM(n, r), and fl25ll26p be- 
tween ASM(7T,, r) and CSM(n, r), the corner sum matrix C which corresponds to the 
edge matrix pair (iJ, V) is given by 



Cij = Hi'j = Viji , for each i,j e [0,n], 



i'=i 



j'=i 



and inversely, 



Hij = Cij — Ci-ij , for each i E [n], j E [0,n] 
Vij = Cij — Cij-i , for each i E [0,n], j E [n]. 



(27) 



(2J 



The set CSM(n, 1) of corner sum matrices for standard alternating sign matrices 
was introduced in IMI, and is also considered in 



and (HI]) 



The corner sum matrix which corresponds to the running example of 
is 

/O 0\ 

1 2 2 2 

112 4 4 

1 2 4 4 6 

2 3 5 6 8 

\0 2 4 6 8 10/ 



(29) 
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Proceeding now to sets of monotone triangles, for n G P and r G N, let MT(?t,, r) be 
the set of all triangular arrays M of the form 



Mil • • • Mir 



M21 



M2,2r 



Mnl 



n,nr 



such that 



• Each entry of M is in [n]. 

• In each row of M, any integer of [n] appears at most r times. 

• Mij < Mi^j+i for each i E [n], j G [ir — 1]. 

• Mifij < Mij < M^i for each i G [n — 1], j G [ir]. 

It follows that the last row of any monotone triangle in MT(n, r) consists of each 
integer of [n] repeated r times. 

It can be checked that there is a bijection between ASM(n, r) and MT(n, r) in 
which the monotone triangle M which corresponds to the higher spin alternating 
sign matrix A is obtained by first using ([8]) to find the vertical edge matrix V which 
corresponds to A, and then placing the integer j Vij times in row i of M, for each 
i,j G [n], with these integers being placed in weakly increasing order along each 
row. (Note that there is alternative bijection in which the horizontal edge matrix H 
which corresponds to A is obtained, and the integer i is then placed Hij times in 
row j of M, for each i,j G [n].) For the inverse mapping, for each i G [0,?7,] and 
j G [n], Vij is set to be the number of times that j occurs in row i of M, and A is 
then obtained from V using 

The set MT{n, 1) of monotone triangles for standard alternating sign matrices was 
introduced in [50], and is also studied in, for example, [33l [Ml [51], IM[ 169]. 

The monotone triangle which corresponds to the running example of and fllip is 



2 3 
13 4 4 
1 2 3 3 5 5 
1 1 2 3 3 4 5 5 
112233445 5. 



(30) 
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Proceeding finally to sets of complementary edge matrix pairs, for n G P and r G N 
we define 



CEM(n,r) 



( ( HiQ . . . Hin\ 



fVoi . . . VonW 



\\-f^nO • • • HnnJ \Vnl ■ ■ ■ VnnJ / 

H2k-lfl = Hn-2k+2,n = 0, l^0,2fc-l = Vn,n-2k+2 = r, for all G [[| 

H2k,0 = Hn-2k+l,n = r, Vo^2k = Vn,n-2k+l = 0, for all G [[f J] 

V,i^ij + Hi^j_i + Vij + Hij = 2r, for all i,j G [n] 



(31) 



It can be seen that there is a bijection between EM(n, r) (and hence ASM(n, r)) and 
CEM{n,r) in which the complementary edge matrix pair {H, V) which corresponds 
to the edge matrix pair [H, V) is given by 



V 



'Hij , i+j odd 
r — Hij , i+j even 

'r-Vij , i+j odd 
Vij , i+j even 



for each z G [n], j G [0,n] 
for each i G [0,n], j G [n]. 



(32) 



The complementary edge matrix pair which corresponds to the running example 
of (BD and ([n]) is 



( 



2 

2 



1 2 

2 



2 
1 

2 10 
1110 2 
2 112 



/2 2 

1 1 

1 1 
1 1 2 



0/ 



2\\ 

2 

2 2 

2 2 



10 10 
V2 2 2/ 



(33) 



In analogy with the association of an edge matrix pair to a configuration of a sta- 
tistical mechanical model, each entry of a complementary edge matrix pair can be 
assigned to an edge of the lattice, i.e., Hij is assigned to the horizontal edge between 
(i, j) and (i, j+l), for each i G [n], j G [0,n], and Vij is assigned to the vertical edge 
between (i,j) and (i + l,j), for each i G [0,n], j G [n]. Also, in analogy with (fT2l) . 
we define the set of complementary vertex types as 



V(r) := {(/i, t;, h\ v) G [0, r]^ | h+v + h'+v = 2r}, 



(34) 
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so that the lattice point is associated with the complementary vertex type 

{Hij_i,Vij, Hij,Vi-ij) G V(r), for each i,j G [n]. Note that the mappings of each 
{h,v, h' ,v') G V(r) to {h,v,r — h' ,r — v'), or of each {h,v, h' ,v') G V(r) to (r — 
h,r — V, h' ,v'), give two bijections between V(r) and V(r). The assignment of the 
entries of the complementary edge matrix pair of ( l33l) to lattice edges is shown 
diagrammatically in Figure [71 
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Figure 7: Assignment of entries of (133!) to lattice edges. 

It is now natural to define, for each n G P and r G N, the set FPL(n, r) of fully 
packed loop configurations as the set of all sets P of nondirected open and closed 
lattice paths such that 

• Successive points on each path of P differ by (—1, 0), (1, 0), (0, —1) or (0, 1). 

• Each edge occupied by a path of P is a horizontal edge between and 

+ with i G [0,n] and j G [n], or a vertical edge between and 
(z + l,j) with i G [n] and j G [0,?7,]. 

• Any two edges occupied successively by a path of P are different. 

• Each edge is occupied by at most r segments of paths of P. 

• Each path of P does not cross itself or any other path of P. 

• Exactly r segments of paths of P pass through each (internal) point of [n]x[n]. 

• At each (external) point (0, 2k — 1) and (ra + l, n — 2k + 2) ioi k e [[f ]], and 
(2/c, 0) and {n — 2k + l,n + l) for k G [[fj], there are exactly r endpoints of 
paths of P, these being the only lattice points which are path endpoints. 

Note that an open nondirected lattice path is a sequence (pi, . . . ,Pm) of points of Z^, 
for some m G P, where the reverse sequence {pm, ■ ■ ■ ,Pi) is regarded as the same 
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path. The endpoints of such a path are pi and pm, and the pairs of successive 
points are pi and p^+i, for each i E [m—l]. A closed nondirected lattice path is a 
sequence (pi, . . . ,Pm) of points of Z^, where reversal and all cyclic permutations of 
the sequence are regarded as the same path. Such a path has no endpoints, and its 
pairs of successive points are pi and p^, for each i G [m — 1], as well as pi and pm- 
For the case of P G FPL(n, r), a path of P whose points are all internal, i.e., in 
[n] X [n] , is closed, and a path of P which has two external points, necessarily its 
endpoints, is open, even if the two external points are the same. 

It can now be seen that there is a mapping from FPL(?2, r) to CEM(n, r) in which 
the fully packed loop configuration P is mapped to the complementary edge matrix 
pair {H, V) according to 

Hij = number of segments of paths of P which occupy the edge between 
(i, j) and + for each z G [n], j G [0,n] 

(35) 

Vij = number of segments of paths of P which occupy the edge between 
+ and for each i G [0,n], j G [n]. 

A fully packed loop configuration of FPL(5, 2) which maps to the complementary 
edge matrix pair of ([3SD is shown diagrammatically in Figure [Hi 
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Figure 8: A fully packed loop configuration which maps to (!33|) . 

It can be checked that the mapping of ( l35l) is surjective for each r G N and n G P. 
Furthermore, for r G {0, 1} or n G {1, 2} it is injective, while for r > 2 and n > 3 
it is not injective. This is due to the fact that if, for a complementary vertex type 
{Hij_i,Vij, Hij,Vi-ij) G V(r), fl35l) is used to assign appropriate numbers of path 
segments to the four edges surrounding the point (i, j) G [n]x[n], then for r G {0, 1} 
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there is always a unique way to link these segments through whereas for 

r > 2 there can be several ways of linking these segments through such cases 

occurring for each n > 3. For example, for r = 2 there is a unique way of linking the 
segments, except if (-^ij-i, V^, ^j-ij, Hy) = (1,1, 1,1), in which case either of the 



configurations 



— or — I — 



can be used. Thus, since the example {H, V) of fl55]) 



and Figure [7] has the single case = (2,2) where this occurs, there are two fully 
packed loop configurations of FPL(5, 2) which map to {H, V): that of Figure [8] and 
that which differs from it by the configuration at (2, 2). 

The cases of FPL(n, 1), and of certain related sets which arise by imposing additional 
symmetry conditions, have been studied extensively. See for example [HI [20l [271 
[28| Wf\ [71]. In these studies, the fully packed loop configurations P of FPL(ra, 1) 
are usually classified according to the link pattern formed among the external points 
by the open paths of P. This then leads to important results and conjectures, 
including unexpected connections with certain statistical mechanical models. See 
for example [23l [23] and references therein. 



Link patterns for certain higher spin statistical mechanical models have been studied 
in [73], although an additional condition is imposed there, which for fully packed loop 
configurations would be that the two endpoints of any open path of P G FPL(?2, r) 
must be distinct. However, the subset of FPL(n, r) in which this condition is satisfied 
is still not in bijection with CEM(n, r), for n > 3 and r > 2 . For example, the 



fully packed loop configurations = 





I 




1 




















1 1 


1 





— and = 



L 



J 



1 



— in FPL(3,2) 



both satisfy the condition and map to the same complementary edge matrix pair in 
CEM(3,2) using ([5^. 



6. The Alternating Sign Matrix Polytope 

2 

In this section, we define the alternating sign matrix polytope in M" , using a half- 
space description, and we show that its vertices are the standard alternating sign 
matrices of size n. 

We begin by summarizing the facts about convex polytopes which will be needed 
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here. For further information, see for example [72j. For m G P, a convex polytope 
in can be defined as a bounded intersection of finitely-many closed affine half- 
spaces in M™, or equivalently as a convex hull of finitely-many points in M™. The 
equivalence of these descriptions is nontrivial and is proved, for example, in [721 
Theorem 1.1]. It follows that hyperplanes in can be included together with 
closed halfspaces in the first description, since a hyperplane is simply the intersec- 
tion of the two closed halfspaces which meet at the hyperplane. The dimension, 
dimP, of a convex polytope V is defined to be the dimension of its affine hull, 
aff(P) := {Api + (1 — A)p2 I Pi,P2 G P, A G M}. A face of V is an intersection 
of V with any hyperplane for which P is a subset of one of the two closed halfs- 
paces determined by the hyperplane. If a face contains only one point, that point 
is known as a vertex. Thus, the set of vertices, verfP, of P C is the set of 
points p & V for which there exists a closed affine halfspace S in M*" such that 
V n S = {p}. It can be shown that verfP is also the set of points p & V which 
do not lie in the interior of any line segment in P, i.e., p G P is a vertex of V if 
and only if there do not exist A G M(o,i) and pi ^ P2 ^ V with p = \pi + (1 — A)p2- 
Any convex polytope V has only finitely-many vertices, and is the convex hull of 
these vertices, or equivalently the set of all convex combinations of its vertices, 
^ = {EpGvertT' -^p P I K ^ ^^[0,1] for cach p G veitV, Y.p&^cvtv\ = !}• Convex 
polytopes V C and V C are defined to be affinely isomorphic if there 
is an affine map (j) : M™' which is bijective between V and V . In such 

cases, verfP' = ^(vert'P). Finally, a convex polytope whose vertices all have integer 
coordinates is known as integral. 

We now define, for n G P, 



■n 



Y.f=iXij' = J21'=iXi'j = 1 for all i,j G [nY 



< 



x = 



fxu ... 




G 



nxn 



J2f=i Xij' > for all i e [n], j e [n — 1] 

YTj'=j > for all i e[n], j e [2, n] 

Y.i'=iXi'j > for all i E [n — 1], j e [n] 

YJi'=iXi'j > for all i G [2,n], j G [n] 



> 



(36) 



= < 



x = 



fxu ... Xin\ 




G W' 



YJj'=iXij' = T,7'=i^i'j = 1 all i,j G [n] ^ 
< J2f=i Xij' < 1 for all i e [n], j e [n-1] 
< Ei'=i Xi'j < 1 for all i e [n-1], j G [n] _ 



21 



In other words, An is the set of n x n real-entry matrices for which all complete row 
and column sums are 1, and all partial row and column sums extending from each 
end of the row or column are nonnegative. 

It can be seen that each entry of any matrix of An is between —1 and 1, and that 
if the entry is in the first or last row or column, then it is between and 1, so 
that An is a bounded subset of M."^ . Since An is also an intersection of finitely-many 

2 2 

closed halfspaces and hyperplanes in , it is a convex polytope in , and will 
be referred to as the alternating sign matrix polytope. This polytope was defined 
independently, using a convex hull description, in [M] . 

An example of an element of A4 is 



/.3 





.6 


.1 


\ 


.2 


.5 


-.6 


.9 


.5 


-.5 


1 







\o 


1 








/ 



(37) 



Defining 

Bn '■= {x G An I Xij > for each i,j G [n]} , (38) 

it can be seen that this is the set of doubly stochastic matrices of size n, i.e., non- 
negative real-entry n x n matrices for which all complete row and column sums 
are 1. This is the convex polytope in M"^ often known as the Birkhoff polytope (or 
assignment polytope). See for example [721 Ex. 0.12] and references therein. 

It now follows that the higher spin alternating sign matrices and semimagic squares 
of size n with line sum r are the integer points of the r-th dilates of An and Bn 
respectively, 

ASM(n, r) = rAn n Z"' , SMS(n, r) = rBn n N"', (39) 

where the r-th dilate of a set P C M*" is simply rP := {rx \ x G P}. 

It also follows that aSAn = aS.Bn = {x e M"^*" | Y.}=iXif = E7'=iXi'j = 1 for all 
i,j G [n]}, and that of the 2n linear equations in variables within this set, only 
2n — l equations are independent, so that 

dimAn = dimBn = {n-lf. (40) 

This is effectively equivalent to the fact that any x G afF^„ = affi3„ can be obtained 
by freely choosing the entries of any (n — l) x (n — l) submatrix of x, the remaining 
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2n—l entries of x then being determined by the condition that each row and column 
sum is 1. 



We also define, for n G P, 

f fhio ... hin\ /f 01 • • • VQn\\ 



ih,v) 



\hnO ■ ■ ■ Kn) \Vnl ■ ■ ■ VnnJ J 



(41) 



hio = voj = 0, h 



1, hij^i+Vij = v,i^ij + hij, for all i,j G [n]\. 



This is a convex polytope in M2n(nrf-i)^ which we shall refer to as the edge matrix 
polytope. It can be seen that EM(n,r) = r£n H Z^"*^"^^\ and that, analogously 
to ([8]) and Qj, there is a bijection between An and £n in which the {h, v) G Sn which 
corresponds to a; G An is given by 



hij = ^ Xij' , for each i G [n], j G [0, n] 



j'=i 



Xj/j , for each i E [0,n], j E [n] 



i'=i 



and inversely. 



Xi 



hij — hij^i = Vij — Vi-ij , for each i,j G [n]. 



(42) 



(43) 



Furthermore, fH2l) can be extended to a linear map from M"' to M2n(r^i)^ g3j) 
can be extended to a linear map from M^'^^"^^) to M"^, implying that An and are 
affinely isomorphic. Also, in analogy with (ITU]) . 



n{n + l)/2 for each (/i,!!) G 



(44) 



It is shown in p!Ol [66] that the vertices of the Birkhoff polytope Bn are the permu- 
tation matrices of size n, so that Bn is an integral convex polytope. We now state 
and prove the corresponding result for An- This result was obtained independently 
in il. 



Theorem 1. The vertices of the alternating sign matrix polytope An are the stan- 
dard alternating sign matrices of size n. 
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Proof. We shall show that the vertices of the edge matrix polytope Sn are the 
edge matrix pairs of EM(n, 1), i.e., vert£^ri = EM(n, 1). It then follows, since Sn 
and An are affinely isomorphic with mapping fH3l) . and since fH3|) maps EM{n, 1) to 
ASM(r;,, 1), that vert^„ = ASM(?t,, 1) as required. 

We first show that EM{n, 1) C vert£:„. Consider any {H, V) G EM{n, 1). From ©, 
this is a pair of matrices with a total of 2?t,(?t, + 1) {0, l}-entries, of which, due 
to ( ITOl) . n{n+l) are O's and n(n+l) are I's. Now define the halfspace S = {{y, z) G 
^nx{rH-i) xM{'>+-i)x" | J2lj=i{Hij Vij + Vij z,j) > ^(^ + 1)}, and consider any matrix 
pair {h,v) G fl 5. Due to ([7]), (14T]) and (jH]), one such matrix pair is {H,V). 
Also, {h,v) G £n implies, using ( HT]) and (jH]), that each of the 2n{n + l) entries 
of {h,v) is between and 1 inclusive, and that they all sum to n{n + l), while 
{h,v) G S implies that the n(n+l) entries of {h,v) in the same positions as the I's 
of {H,V) sum to at least n{n+l). It can be seen that these conditions are only 
satisfied if {h,v) = {H, V). Therefore, SnHS = {{H, V)}, implying that {H, V) G 
vert^n as required. (Note that alternatively it could have been shown here that 
G M"x(,^i)x^(rH-i)xn | YJl^.^i Hij > n{n + l)/2} = {iH,V)} or that 

Sn n {{y, Z) G M"X(rH-l) xM(-l-l)xn | ^n^^^ y,^ > n{n + l)/2} = {{H, V)}.) 

We now show that vert£^„ C EM(n, 1). Consider any {h, v) G £^ri\EM(n, 1). We shall 
eventually deduce that {h,v) ^ vert£^n, which gives the required result. Similarly to 
the association of edge matrix pairs with configurations of a statistical mechanical 
model, we associate with the horizontal edge between lattice points and 
+ for each i G [n], j G [0,n], and Vij with the vertical edge between lattice 
points and + for each i G [0,n], j G [n] (using matrix-type labeling of 
lattice points). Since EM(n, 1) = ^„ n Z2"('^i), {h,v) ^ EM(n, 1) implies that at 
least one entry of {h,v) is nonintegral. Now, from fHTl) . {h,v) G Sn implies that 
hio = Vqj = and hin = Vnj = 1, for each i,j G [n], so that any nonintegral entry of 
{h, v) must be associated with one of the 2r;,(r;,— 1) internal edges, (i.e., the horizontal 
edges between and (i, j+1), for i E [n], j E [n—1], and the vertical edges between 
and + for each i G [^ — 1], j G [n]). Also from (HTl) . {h,v) G Sn implies, 

that 

hij_i + Vij = Vi^ij + hij, (45) 

for each i,j G [n]. But if any one of the four entries in fl45p is nonintegral, then at 
least one of the others must also be nonintegral. Therefore, the existence among 
the entries of {h, v) of a noninteger implies the existence of two further nonintegers, 
among each of the other three entries of the two cases of fH5|) in which the initial 
nonintegral entry appears. It now follows by repeatedly applying this argument. 
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and since the internal edges form a finite and closed grid, that there exists at least 
one cycle of internal edges associated with noninteger entries of {h,v). 



We select any such cycle, give it an orientation, say anticlockwise, and denote the 
sets of points (i, j) for which the horizontal edge between and is in the 

cycle and directed right or left as respectively 7i+ or 7i_, and the sets of points 
for which the vertical edge between and + is in the cycle and directed 
up or down as respectively V+ or An example of such a cycle, for the {h, v) G £4, 
which corresponds to the example x G ^4 of fl57|) is shown diagrammatically in 
Figured For this example, n+ = {(2, 2), (2, 3), (3, 1)}, = {(1, 1), (1, 2), (1, 3)}, 
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Figure 9: A cycle of nonintegers for the example of ( 137|) . 



V+ = {(1,4), (2,2)} and V_ = {(1, 1), (2, 1)}. We now define, for ^ G M, a matrix 
pair (/i'(/i), t>'(yu)), with entries h'{^)ij for i G [n], j G [0,n], and v\ji)ij for i G [0,n], 
i G [n], given by 



hij+fi, (i, j) G n+ 
/i, (ij) G H- 
otherwise 



hij ■ 
h 



Vij-^^, (ij) G V_ 
otherwise . 



For the example of Figure M, 
I 

{h{fi),v{fi)) = 



/O 




Vo 



.3 — /i 
.2 





.3 — /i 
.7+/i 



1 



.9-/i 

.l+^J. 
1 
1 



1 
1 

V 



/ 

3 — /X 
5 — ^ 
1 

V 1 






.5+1-1 

1 




.6 


1 
1 



\\ 

+ 
1 
1 

1 // 



(46) 



(47) 



We now check whether {h' (fi) , v' (fi)) G Sn- By using (l46l) to replace each entry of 
{h,v) in psl) with an entry of (/i'(/i), f '(/i)), it can be checked that the equation 
h'{fi)ij^i+v'{fi)ij = v'{fi)i^ij+h'{fi)ij is satisfied for each i,j G [n], since if the cycle 
does not pass through (i, j), then the required equation is immediately obtained. 
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while for all possible configurations (of which there are six), and both possible 
directions, in which the cycle can pass through all explicit appearances of /i 

cancel out, again leaving the required equation. The conditions h'{fi)iQ = t>'(/i)oj = 
and h'{fi)in = v'{fi)nj = 1 are also met for each i,j G [n], since (/i'(yu), and 
{h,v) match for these entries. It only remains for the conditions h'{fi)ij,v'{^)ij G 
]R[o,i] to be checked for each i,j G [n], and it can be seen that these are satisfied if 
— /i„ < ^ < ^1+1 where 

/i_ := min({/iij | (i, j) G 71+} U {vij \ G V+} 

U {1-hij I G H-} U {l-v,j I G V-}) 

(48) 

/i+ := mm{{l-hij \ G n+} U {!-%• | G V+} 
U {hj I G U I (z, j) G V-}) . 

The facts that /ijj G K(o,i) for each (i, j) G TY-t, and fjj G IR(o,i) for each (z, j) G V±, 
imply that yU_ and yU+ are both positive, so that there exists a finite-length, closed 
interval of values of /i for which {h' (fi) , v' (fi)) G Sn- For the example of Figure [HI 
/i_ = .1 and /i+ = .3. Finally, we define {h±,v±) := (/i'(±/i±), t;'(±/i±)), so that it 
follows that {h±,v±) G £n, {h_,vJ) ^ {h+,v+), and 

(h, ^) = ^(^-' + ^(^+' • (49) 
Therefore, {h, v) lies in the interior of a line segment between two points of Sn, and 
so, as required, {h,v) ^ vert£^„. This concludes the proof of Theorem 1. □ 

It follows immediately from Theorem 1 that the alternating sign matrix polytope is 
integral, a fact which will be important in the enumeration of higher spin alternating 
sign matrices in the next section. 

It also follows from Theorem 1 that An can be described as the convex hull of the 
alternating sign matrices of size n, or equivalently that 

A = 

\a G M[o,i] for each A G ASM(n, 1), ^ = 1 1. 

AGASM(n,l) J 

The way in which this conclusion has been reached here depends on the general 
theorem that any convex polytope has both a halfspace and convex hull description. 
However, fl50p can also be derived more directly. It follows immediately from ([T]) 
and fl36|) that the RHS of fl50l) is a subset of the LHS, i.e., that every convex com- 
bination of standard alternating sign matrices of size n is an element of An- Con- 
versely, every element of An is a convex combination of elements of ASM(n, 1) due 



lAeASM(n,l) 
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to the following argument, which we owe to [M]. For each x G An, let z/(a;) be the 
number of nonintegers among the entries of {h,v), where {h,v) G Sn corresponds 
to X by P2|) . Now consider the following recursive process applied to x G An- If 
z/(a;) = 0, then x G ASM(n, 1), and the process is terminated. If z/(a;) is positive, 
then let x± G An correspond to (/i±, f ±) G £n, where {h±,v±) are defined in terms of 
{h,v) as in the proof of Theorem 1. Using (l49l) . write x as the convex combination 

~ JT^+JI^-'^- /7"^Tm+"^+' ^^^^^ ^^^^ ^^^^ '^(^^i) < '^(a;)- Now reapply 

the process to x_ and x^. It follows that any x G ^„ can thus be completely 
decomposed as a convex combination of elements of ASM(n, 1). 

Faces of the alternating sign polytope other than those given by its vertices are 
studied in detail in [M], but we shall not consider these here. 



7. Enumeration of Higher Spin Alternating Sign Matrices of Fixed Size 

In this section, we use the general theory of the enumeration of integer points in 
integer dilates of integral convex polytopes to obtain results on the enumeration of 
higher spin alternating sign matrices of fixed size. 

We begin by summarizing the details of this theory which will be needed here. For 
further information, more general results, and references, see for example [7] or [62l 
Sec. 4.6]. For an integral convex polytope V in with relative interior V°, there 
exists a unique function Cp{r), where r can be regarded as a complex variable, with 
the properties that: 

(i) C-plr) is a polynomial in r of degree dim P. 



(ii) \rV n Z™| = Cvir), for each r G N. 
^iii) \rV° n Z"^| = (-1)'^''^^ ^^(-r), for each r G N 



(iv) C-pir) can be expressed in the form C-pi^r) 



(51) 



dimP / 



dimP 



with Cfc G N for each G [0, dimP]. 

The function Cp{r) is known as the Ehrhart polynomial of V. It can be seen that 
£p(0) = Cdim-p = 1, and that the Ehrhart polynomial can be obtained explicitly by 
finding >C-p(r) for dimP further integers r by directly enumerating the number of 
lattice points \rV fl Z*"! for positive r, or | — rV° fl Z^l for negative r, and then 
interpolating. 
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For the alternating sign matrix polytope An, the relative interior A° is obtained by 
simply replacing each weak inequality in (|36|) by a strict inequality. Defining, for 
n G P and r G N, 

ASM° (n,r) := nZ"', 



it is seen that 
ASM°(n,r) := 



Mn ... AiA 

\Anl ■ ■ ■ Ann/ 



E/=iA^/ = J:7'=iA'j = r for all t,j G [n]] 

v4jj/ > 1 for all i e[n], j e [n-1] 
Ef=j Af > 1 for all 2 G M, J G [2, n] 
Ei'=i A'j > 1 for all i G [n-1], j G [n] 
E7'=iAi'j > 1 for all i G [2,n], J G [n] 



(52) 



(53) 



Mil . 


■ Am^ 




• Ej'=iAij' — T,7'=iA'j 


= r for all i 


J e [n] 


-,i 






G Z"^" 


• 1 < Ef=iAiy<r-l 


for all z G [n 


, j e [ri- 






■ Ann J 




•l<j:i'=iAi'j <r-l 


for all i G [n- 


-1], je 





Thus, ASM° (n,r) is the set of nxra integer-entry matrices for which all complete 
row and column sums are r, and all partial row and column sums extending from 
each end of the row or column are positive. It follows that 



A^ni + l 
and that 



Aln + l\ 
Ann + 1/ 



/Au ...Ain\ 

\Anl ■ ■ ■ Ann) 



G ASM(n,r-n) 



\ C ASM°(n,r) 



(54) 



for each n, r G P with r > n, 



ASM°(n,r) 



for each n, r G P with r < n. 



(55) 



It can be seen immediately that the containment of flMj) is in fact an equality for 
n = 1 and n = 2, and it follows from flBBl) in Section 8, that this is also the case for 
n = 3. 

We now state and prove the main result for the enumeration of ASM(n, r) and 
ASM°(n,r). 
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Theorem 2. For fixed n G P, there exists a function A„(r), the Ehrhart polynomial 
of the alternating sign matrix polytope An, which satisfies: 

(i) An{r) is a polynomial in r of degree (n — 1)^. 
(a) |ASM(n,r)| = v4„(r), for eachr eN. 
(ill) |ASM°(n,r)| = (-l)"+M„(-r), for eachr en. 
(iv) An{~l) = Ani-2) = ... =Ani-n+l) = 0. 

(vi) An{r) can he expressed in the form A„(r) = 

fc=n-l \ 



\n-l) 



2 



wi/i Cfc G N /or each k G [ri — 1, (ri — 1) 



Proof. All of the conclusions of this theorem follow straightforwardly from results 
already obtained or stated in this paper. First, the existence of an Ehrhart poly- 
nomial Cj[^{r) = An{r) follows from the fact, implied by Theorem 1, that An is 
an integral convex polytope. Properties (i)-(iii) in this theorem then follow from 
(i)-(iii) of (15T]) applied to An, using (139|) . (l40l) and (!52|) . Property (iv) follows from 
property (iii) applied to (1551) . Property (v) follows from property (ii) applied to (I6l). 
Finally, property (vi) follows from (iv) of 0511) . and from property (iv), which gives 

Co = Ci = . . . = Cn-2 = 0. □ 

It follows that the explicit polynomial A„(r) for a particular n G P can be found 
by interpolation using the n + 1 values provided by An{0) = 1 and properties (iv) 
and (v), together with n^— 3n+l further values obtained by the direct enumeration 
of cases of ASM(?T,,r) or ASM°(?T,,r) and the application of properties (ii) and (iii). 
We have done this for n = 3, 4 and 5, some of the required values being provided in 
Table [H Together with the trivial cases n = 1 and 2, the results, expressed in the 
form of (vi) of Theorem 2, are 



9 / V 9 



(57) 
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and 



Adr) = 70^+/) + 14468(" + ^) + 52165l(" + ^) + 6002192(" + ^^ + 
282335651 ' ] + 61083124 ( _ + 64066830' 



16 y V 16 



32866092 f^ti^l + 7998192^13^^ +854464^^3^') + 



16 y V 16 / V 16 

^r + 14\ ,,„/r + 15\ /r+16~ 



16 



(5J 



We note that the fact that fl5^ is an equahty for n = 1, 2 and 3, imphes that 
|ASM°(n, r)| = |ASM(?T,, r— for all n G {1, 2, 3} and r G P with r > n, which in 
turn implies that A„(r) = {—l)"'^^An{—n—r), for all n G {1, 2, 3} and r G C, as can 
be observed in fISBl) . 

Finally, we now outline the results for Bn and SMS(n, r) which correspond to those of 
this section for An and ASM(n, r). Even though the results for Bn and SMS(ra, r) are 
previously-known, having been conjectured in [2] and first proved in [311 EI], we state 
them here to show their similarity with those for An and ASM(r;,, r). The relative 
interior 0° of the Birkhoff polytope Bn is the set of positive real-entry nxn matrices 
for which all complete row and column sums are 1, and SMS°(n, r) := ri3° fl Z"^ is 
the set of positive integer-entry matrices for which all complete row and column sums 
are r. It follows that if SMS(n, r) and SMS°(n, r) are substituted for ASM(n, r) and 
ASM°(n, r) in fl54p and fl55p . then the equations still hold, and furthermore that the 
containment in fl5^ can be replaced by an equality in all cases. Since Bn is an integral 
convex polytope, it has an Ehrhart polynomial Ci3„{r), which we denote Hn{r). 
Applying flSTl) . and using ([5]), fl39l) and fHOl) . it now follows that all of Theorem 2 still 
holds if the substitutions Hn{r) for v4„(r), 'Birkhoff polytope' for 'alternating sign 
matrix polytope', Bn for X, SMS(n,r) for ASM(n,r), SMS°(n,r) for ASM°(n,r), 
and n\ for nr=o(3^+l)'/('^+^)' made. Furthermore, the fact that the counterpart 
to (154]) is now an equality implies that |SMS°(n, r)| = |SMS(n, r— for all n, r G P 
with r > n, which leads to the further properties that Hn{r) = (— l)"+^if„(— n— r) for 
all n G P and r G C, and, in the counterpart to (vi) of Theorem 2, that = Cn{n-i)-k 
for each n G P and k E [n — 1, (n — 1)^]. 
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8. Higher Spin Alternating Sign Matrices of Size 3 



In the previous section, the enumeration of higher spin alternating sign matrices 
was studied using a general, but nondirect, approach. In this section, we consider 
the special case of 3 x 3 higher spin alternating sign matrices, and provide a direct 
bijective derivation of the enumeration formula, from Theorem 2 and fl56|) . 

|ASM(3,r)| = (^^'^^ + ("^t^' '''^''^ ^ ^' ^^^^ 



We begin by outlining a derivation of the corresponding semimagic squares formula, 

|SMS(3,r)| = (^^ + 2) + + + ('■+^), for each r G N, (60) 

since the derivation of (l59l) will be similar. The polynomial fl60|) for |SMS(3,r)| was 
obtained directly in [H [21 [121 [ISl [S2] • The derivation given here most closely follows 
that of 



Let 5*3 be the set of permutations of 123, and consider the set 
C{r) : = 

• J2aeS3 O'er = r 



(,Ol23; ^132, 0213; 

0231, «312, 0321) G 



(61) 

ai32 = or (a32i = and 0132 7^ 0) ^ • 
or (a2i3 = and 0132, 0321 7^ 0) 
Using the fact that, for any m, n G N, the number of tuples (oi, . . . , cim+n) G P^xN" 

with y]™it" ttj = r is ^ ) (which can be obtained by considering tuples 

\m+n — lj 

(61, ... , bm+n~i) e P™+"-i with bi <...<bm< bm+i <■■■< bm+n-i < r, and using 
the bijection 6j = I]}=i (^j for each i G [m+n—1]), it follows that 

icwi^rrurt^um. (62) 



. 4 y V 4 y V 4 . 

where the three binomial coefficients correspond respectively to the three alterna- 
tives in the second condition of (1611) . 

Now let Pa be the permutation matrix given by {Pa)ij '■= Si,aj for each cr G 6*3 and 
i,j G [3], and define, for each a G C(r), 

/ctl23 + Ctl32 ^213 + 0312 a231 + '^321^ 



9ia) := ^ a^Pa 



0213 + 0231 0123 + 0321 O132 + O312 I . (63) 
V'^312 + 0321 O132 + O231 Oi23 + 0213y 
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It can immediately be seen that 6{a) G SMS(3,r) for each a G C{r). It can also 
be checked straightforwardly that 6* is a bijection between C{r) and SMS(3,r), the 
inverse mapping being, for each A G SMS(3,r), 



■ (An, 0, ^33-^11, A32, A23, ^22-^11), 

e~\A) = I (^22,^11-^22,^33-^2,^13, ^31,0), 



(A33, An - A33, 0, A21, A12, A22-A 



33 J 



All < A22 and All < A33 
A22 < All and A22 < A33 
A33 < All and A33 < A22, 



(64) 

where the three cases of (1641) correspond respectively to the three alternatives in the 
second condition of (!6T|) . Therefore |SMS(3,r)| = |C(r)|, so that (!62|) gives (1601) . 
thereby completing its bijective derivation. 

Proceeding to 3x3 higher spin alternating sign matrices, we first note that ASM(3, r) 
is simply the set of 3 x 3 integer-entry matrices for which all complete row and column 
sums are r, and the eight external entries are nonnegative. 



ASM(3,r) 





(All 


A12 


Al3\ 


• 




A21 


A22 


A23 G 






[A31 


A32 


A33J 


• 



E-'=i^y' = Ei'=iAi'j = r for all z, j G [3] | (65) 
A,, >0 for all (^,j) e [3]x[3]/{(2,2)} 



To obtain fl65|) from ([T]), it only needs to be verified that if A is an element of the 
RHS of ([65]), then each partial sum A12 + ^22, A21 + ^22, A23 + A22 and ^132 + A22 
is nonnegative. For the case of A12 + A22, ^12 + ^22 + ^32 = 7431 + 7432 + ^33 (= r) 
gives A12 + A22 = 7431 + ^433, but A31 > and ^33 > now imply that A12 + A22 > 
as required. The other three cases follow by symmetry. 

The set of standard alternating sign matrices, ASM(3, 1), is the set {P^ \ a G S3} 



of 3 X 3 permutation matrices, together with the matrix 




which will be 



denoted as Pq. Setting 
C'{r) := 



I ._ 



5*3 U {0}, we now define, in analogy with (IMl) and (1^ . 



(ai23, ai32, ^213, (^231, 

0312,0321,00) G N'^ 



Eaes'^ aa = r 

Oi32 = oo = or 

(0321 = oo = and 0132 7^ 0) or 

(0213 = ao = and 0132, 0321 7^ 0) or 

(0123 = 0321 = and ao 7^ 0) 



(66) 
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0.123 + 0132 O213 + O312 + O0 0231+0321 
0213 + 0231 + ^0 Oi23 + O321~0o 0132 + 0312 + 00 
O312 + O321 O132 + O231 + O0 O123 + O213 



(67) 



cv)i^CtVKT)^(T) 



(68) 



and that 9'{a) G ASM(3,r) for each a G C"(r). It can also be checked, using (1^ 
and the bijection fl63|) and fl64l) between C(r) and SMS(3,r), that 9' is a bijection 
between C'{r) and ASM(3,r), the inverse mapping being, for each A G ASM(3,r), 



'(All, 0,^33-^11,^32,^23,^22-^11,0), Al < ^22 and Au < A33 

(^22,^11-^22,^33-^22,^3, ^31, 0,0), < A22 < and A22 < A33 (69) 



(^33,^11-^33, 0,^21, ^12,^22-^33,0), A33 < Au and A33 < A22 

^ (0, All, A33, Ai3, A31, 0, -A22), A22 < 0, 

where the four cases of fl69l) correspond respectively to the four alternatives in the 
second condition of fl66|) . Therefore |ASM(3,r)| = |C"(r)|, so that fl68l) gives fl59l) . 
thereby completing its bijective derivation. 

In this section, particular sets C{r) and C'{r) of weak compositions of r into six 
and seven parts, with at most five parts nonzero, have been found which are in 
bijection with SMS(3, r) and ASM(3, r) respectively, where the mapping from a weak 
composition to a matrix simply uses the parts of the composition as coefficients in 
a linear combination of permutation matrices or standard alternating sign matrices. 
It would clearly be worth investigating whether these results can be generalized to 
larger integers n to give sets of weak compositions of r into n! and Yi^Zoi^i+^V- / 
parts, with at most (n — 1)^ + 1 parts nonzero, which are similarly in bijection with 
SMS(n,r) and ASM(n, r) respectively. This, however, seems to be a challenging 
problem. 

9. Discussion 

We conclude this paper by outlining some possible directions for further research. 
While all of these can be simply stated, and follow from the content of this paper 



9'-\A) 
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in obvious ways, some may well turn out to be of limited interest, or to lead to 
prohibitively-difficult problems. 

In this paper, we studied the configurations of certain integrable statistical mechan- 
ical models with domain wall boundary conditions. One area for additional work 
would therefore be the consideration of further integrable statistical mechanical mod- 
els with boundary conditions of domain- wall- type. A possible component of such 
work would be the derivation of determinant formulae for the partition functions of 
such models, this having already been done for certain models in [TS| [2^ . 

In Section 3, we discussed the determinant formula for higher spin vertex models. 
Although it seems unlikely that this formula can be used for the unweighted enu- 
meration of higher spin alternating sign matrices or semimagic squares, it may still 
be possible for it to be used to obtain certain weighted enumeration formulae. More 
specifically, each vertex type of V(r), as defined in (fT2|) . would be given a weight, 
and the overall weight of a higher spin alternating sign matrix would be the prod- 
uct of the weights of all the vertex types associated with the corresponding edge 
matrix pair. Certain such weighted enumeration formulae are already known for 
standard alternating sign matrices (see for example [211 E21 US])- In these cases, 
vertex type (4) in Figure [3] (which corresponds to each entry of —1 in a standard 
alternating sign matrix) is given a weight of 2 or 3, and the other five vertex types 
of V(l) are each given a weight of 1. 

In Sections 4 and 5, we considered various other combinatorial objects, namely 
certain sets of lattice paths, monotone triangles, corner sum matrices and comple- 
mentary edge matrix pairs, which are in simple bijection with higher spin alternating 
sign matrices. Each of these representations could be studied further, and might pro- 
vide useful statistics according to which higher spin alternating sign matrices could 
be weighted and classified. The fully packed loop configurations defined in Section 5 
also seem worthy of further study, especially since their classification according to 
the link pattern formed among the external points by the open paths, may provide 
an interesting generalization of such a classification for the standard case of r = 1 
(see for example [ISl EDI EZl [2HI EZl El] ) • 

It is also common to impose further conditions on semimagic squares (see for ex- 
ample [6] and [621 Proposition 4.6.21]) and standard alternating sign matrices (see 
for example [IHl E21 ESI ESI EZl)- It would therefore seem natural to investigate the 
effects of imposing similar such conditions on higher spin alternating sign matrices. 
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it being expected that this would lead to the introduction of further, not necessarily 
integral, convex polytopes. 

Semimagic squares are a special case of contingency or frequency tables (see for 
example [37j), these simply being nonnegative integer-entry rectangular matrices 
with arbitrary prescribed row and column sums. Similar generalizations of higher 
spin alternating sign matrices in which the complete row and column sums are 
prescribed, but not necessarily all equal, and the matrices are not necessarily square, 
would thus provide another direction for further study. 

In Section 6, we defined the alternating sign matrix polytope, and showed that 
its vertices are standard alternating sign matrices. Many additional results for 
this polytope, including the enumeration of its facets, the characterization of its 
face lattice, and the description of its projection to the permutohedron have been 
obtained in [M], but this polytope could still be investigated further. 

In Section 7, we obtained certain enumerative formulae for higher spin alternating 
sign matrices with fixed size n and variable line-sum r, but we have not obtained 
formulae for the case of variable n and fixed r. Formulae for semimagic squares with 
variable n, and r = 2 or r = 3 are known (see for example [2], [22l Ex. 25, Page 124] 
and [63l Sec. 5.5], but note that the r = 3 formula in [22] contains misprints). 
However, it might be expected that the difference in difficulties in obtaining, say, 
the r = 2 semimagic square formula and an r = 2 higher spin alternating sign 
matrix formula, might be comparable to the substantial difference in difficulties in 
obtaining the r = 1 formulae ([5]) and ([6]). 

Finally, as indicated at the end of Section 8, it would be interesting to see whether 
the approach used there can be generalized to give bijective derivations of enumer- 
ative formulae for semimagic squares or higher spin alternating sign matrices with 
fixed size larger than 3. 
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